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The 3d electronic structure and phase transition in pure and Cr doped V2O3 are theoretically 
investigated in relation to the 3d spin-orbit interaction and lattice distortion. A model consisting 
of the nearest-neighbor V ion pair with full degeneracy of the 3d orbitals is studied within the 
many-body point of view. It is shown that each V ion with S — 1 spin state has a large orbital 
magnetic moment ~ 0.7/iB and no orbital ordering occurs in the antiferromagnetic insulating 
(AFI) phase. The anomalous resonant Bragg reflection found in the AFI phase is attributed to 
the magnetic ordering. In the AFI and paramagnetic insulating (PI) phases, Jahn- Teller like 
lattice instability leads to tilting of the V ion pairs from the corundum c-axis and this causes 
large difference in the orbital occupation between the paramagnetic metal and the insulating 
phases, which is consistent with linear dichroic V 2p XAS measurements. To understand the 
AFI to PI transition, a model spin Hamiltonian is also proposed. The transition is found to be 
simultaneous order-disorder transition of the magnetic moments and tilting directions of the V 
ion pairs. Softening of elastic constant C44 and abrupt change in short range spin correlations 
observed at the transition are also explained. 

KEYWORDS: V2O3, electronic structure, phase transition, lattice distortion, orbital ordering, spin-orbit interac- 
tion, orbital moment, resonant x-ray scattering, V 2p XAS 
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§1. Introduction 

The metal-insulator (Mott) transition in V2O3 have 
been extensively studied ever since Mott suggested the 
possibilitjj-pf the locahzation of electrons driven by Coulomb 
repulsion.El Above the transition temperature Tn ~ 
155K, the material is paramagnetic metal (PM) and has 
the corundum structure with the trigonal symmetry (space 
group i?3c). Below Tn, it is an antiferromagnetic in- 
sulator (AFI) and the lattice distortion to the mono- 
clinic structure {12 /a) takes place.© Cr dopecLajkiys 
(Vi_2;Cra;)203 exhibit also interesting features. clQla'B'EP 
For 0.005 < a: < 0.018, at a temperature above the AFI 
— > PM transition, the material further undergoes a PM 
to paramagnetic insulator (PI) phase transition with no 
change in the crystal structure. For much higher Cr con- 
centration X > 0.018, the PM phase disappears and a 
AFI PI transition takes place. 

Although V2O3 and its alloys are often referred to 
a typical example of the Mott-Hubbard system, there 
are experimental facts showing the complex aspects of 
these materials, which can not be explained by the sim- 
ple single-band Hubbard model. In the AFI phase, they 
exhibit unusual mametic order, which breaks corundum 
trigonal symmetry .B^ Neutron scattering experiments 
have been performed by Bao et al. and abrupt change in 
spatial spin correlation functions has been reported hott 
at the AFI PM and the AFI -> PI transitions Ja'EJ 
This switching of spin correlation indicates that the phase 
transitions can not be considered as usual order-disorder 
magnetic transition. Softening of the C44|-elastic con- 
stant has been observed in ultrasonic wavetJ^ and neu- 



tron scatterin^ll^ experiments at the AFI — > PI transi- 
tion, indicating nearly second order feature of this phase 
transition as well as strong electron-lattice coupling. 

Even now, the electronic structure of these materials 
are still controversial problem and above experimental 
facts suggest that much realistic model including the or- 
bital degeneracy of the 3d orbital and the lattice dis- 
tortion is necessary to explain them. In V2O3, the Oh 
crystal field splits the 3d orbitals of each V ion into two 
fold eg and three fold t2g orbitals and the t2g orbitals are 
further separated into low laying two fold e'^ and non- 
degenerate fli orbitals by small trigonal field. In local- 
ization limit, two 3d electrons in the ion are accom- 
modated in the lowest orbitals with = 1 spin state. 
Castellani et al. haiie proposed a model including these 
and ai orbitals £3* In the corundum structure, two 

V ions forms pairs along the c-axis and in their model, 
the oi orbitals on two V ions of the pair forms molecu- 
lar orbitals. The bonding ai orbital accommodates two 
electrons out of four in the pair and remaining two elec- 
trons occupies one of e'^ orbital in each V ions, resulting 
in e'^ai configuration with 5 = 1/2 spin state in each 

V ion. They also propose a long-range ordering of or- 
bital occupation in the two fold orbital to explain the 
complex antiferromagnetic ordering in the AFI phase. 

However, recent measurements on the linear dichroism 
of V 2p x-ray absorption spectra (XAS) in V2O3 have 
revealed that the 3d electron configuration is not S ~ 1/2 
e" ai but mixtucfi of e^e^ and e" ai with 5 = 1 Hund's 
rule spin state.M In addition, large differences in e^eT' 
to e^ai ratio have been found among the AFI, PI and 
PM phases. Indeed, such 5=1 Hund's rule spin state is 
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also obtained in recent "LDA+[/"lla't3 band structure 
calculationJlH-' in which a S* = 1 e^e^ ground state having 
experimentally observed antiferromagnetic order with no 
orbital degeneracy is favorable for the low temperature 
monoclinic phase. Mila et al. have proposed_a, jnodel 
with a large intraatomic exchange interaction Jl3'H In 
this model, spin state in each V ion remains being S — 
1 even with strong hybridization between ai orbitals. 
The electronic state of the V ion pair is described as 
a superposition of two configurations with equal weight: 
e^ai;e'^e'" (the e'^ai configuration in one of V ion and 
the e^e'^ configuration in the other) and e'^e'^; e'^ai. In 
this model, two fold orbital degeneracy is present in this 
V2 molecular orbital and the complex antiferromagnetic 
order in the AFI phase is stabilized with a ferro-type 
order of this V2 molecular orbital. 

Paolasini et al. have performed resonant x-ray scat- 
tering experiments (RXS) at the V K-edge on V2O3 and 
they insisted that resonant Bragg peaks obseived in the 
AFI phase are due to the 3d orbital ordering. EHlElP How- 
ever, Lovesey and Knight have argued that the resonant 
Bragg peaks arise fropLthe ordered orbital magnetic mo- 
ments of the V ionsE3) Indeed, nonresonant magnetic 
scattering has been also measured and a large contri- 
bution from the orbital magnetic moment to tbr; total 
magnetic moment Ml /Ms ~ —0.3 was found.EiJ' It is 
thus important to know how the 3d spin-orbit interaction 
affects the 2>d electronic structure of V2O3. Since this in- 
teraction couples spin and orbital degrees of freedom, in 
the presence of a large orbital angular momentum, no 
orbital ordering is expected in the AFI phase. In above 
mentioned theories on the electronic structure of V2O3, 
however, the effects of the id spin-orbit interaction in 
the V ions is presumed to be small. 

The purpose of this paper is to clarify the effects of the 
3c? spin-orbit interaction and the lattice distortion on the 
Zd electronic structure and the phase transitions of the 
pure and Cr doped V2O3. To this end, first of all, a V 
ion pair model similar to that proposed by Mila et al. 
will be considered taking into account the Zd spin-orbit 
interaction in § 2. With the 3d spin-orbit interaction, 
no orbital ordering is expected in the AFI phase and in- 
stead of this, a large orbital magnetic moment ^ 0.7/iB 
is induced in the ground state. The relation between the 
monoclinic distortion and the magnetic anisotropy and 
difference in the configuration of the 3c? electrons between 
the AFI and PM phases will be also discussed. Second, 
interaction between the electronic state of the V ion pair 
and local lattice distortion caused by displacement of the 
pair will be discussed in § 3. It will be shown that Jahn- 
Teller like lattice instability occurs when the energies 
of the e^ai^e^e" {e'^ e'^ ; e^ ai) and e'^e'^;e'^e^ states are 
nearly degenerate. Third, relation between the magnetic 
ordering and the lattice distortion in the AFI phase and 
also the mechanism of the AFI PI transition will be 
discussed in § 4. For this purpose, an effective Hamilto- 
nian with exchange interaction between the V ion pairs 
including lattice distortion effects is proposed. Softening 
of the elastic constant C44 and abrupt change in short 
range spin correlation functions observed at the AFI — > 
PI transition are reproduced with the model. Finally, to 



confirm the validity of the present model, recent experi- 
ments on RXS and the linear dichroic V 2p XAS spectra 
are analyzed in § 5 and § 6, respectively. The experimen- 
tal results are well explained within the present model. 

§2. V ion Pair Model 

2.1 Model Hamiltonian 
Each V ion in the corundum structure of V2O3 is octa- 




Fig. 1. Corundum structure of V2O3. Each V ion (open cir- 
cles) is octahedrally surrounded by six oxygen ions (octahedra). 
Along the cji-axis, the V ion pairs surrounded by the face-sharing 
two oxygen-octahedra are seen whereas in a^-fth plane, each V 
ion has three neighbors with edge-sharing oxygen-octahedra. 



hedrally surrounded by six oxygen ions (see Fig. 1). An 
octahedral crystal field and hybridization between the 
V 3c? and O 2p orbitals split the energy level of the 3d 
orbital into the eg and t2g levels and the low laying t2g 
orbitals accommodate two electrons in each V'^+ ion. A 
small trigonal distortion of the oxygen octahedron along 
the c^-axis further splits the t2g level into a doubly de- 
generate e'^ and a non-degenerate ai levels with the C^v 
symmetry. In the corundum structure, the layers of the 
honeycomb lattice of V ions are stacking along the Ch- 
axis and each V ion in the layers forms a V ion pair along 
the c;i-axis with a V ion in the adjacent layer. While the 
fli orbitals point to the direction of the cj^-axis, form- 
ing strong a bonds in the V ion pairs, the e'^ orbitals 
extend their lobes toward the direction perpendicular to 
the c^i-axis, having tt bonds with the three neighbors in 
the layer. 

As was first suggested by Allen, here, it is assumed 
that the electrons are delocalized within the nearest- 
neighbor V ion pair mid interaction between much fa- 
ther V ions is week.t3^ The two V ions of the pair 
along the c/j-axis and electron hopping between them 
were coasideced based on configuration interaction ap- 
proachBj'E^Ea' In each V ion site, the fivefold 3d or- 
bitals are fully included and the Coulomb and exchange 
interactions among 3d electrons, the 3d spin-orbit inter- 
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action and the trigonal crystal field are considered. EIP 
The Hamiltonian for the V ion pair model then would 
be written as 

iy\ ,. . . ,i-'4 

+ C^{v,v')al^a,y, (1) 

where the first and second terms are those for the one 
body energies of the 3c? orbitals, the effective Coulomb re- 
pulsion energy between the 3d electrons, respectively, the 
third term represents the electron hopping between two 
V ion sites A and B in the V ion pair given in terms of the 
Slater-Koster parameters dda and ddTr, the forth term is 
for the multipole interaction between the 3d electrons, 
which is described by the Slater integrals Fj^ and F'^^ 
of the 3d orbital and the fifth term denotes the 3d spin- 
orbit interaction with the coupling constant C,d- The last 
term represents crystal field due to the monoclinic lat- 
tice distortion in the AFI phase. Here, v represents the 
combined index of the species of the 3d orbital and the 
spin quantum number of the electron, i = A, B stands 
for the V ion site and aj^^ creates an electron on the 3d 
orbital with v at the V ion site i. 

In a C-iv crystal field, the five degenerate 3d orbitals 
split into one non-degenerate oi orbital and two doubly- 
degenerate e orbitals. Using coordinates where z-axis is 
parallel to the c^-axis of the corundum structure and 
y-axis to one of the three V-V bond direction in the a^- 
bh basal plane (see Fig. 1), the wavefunctions for the 3d 
orbitals in terms of the C^y basis functions are written 
as 



\ai:i) - \3z^~r^), 




where the two kinds of e orbitals in the Csv symmetry 
originating from the t2g (eg) orbitals in the Oh symme- 
try is specified by superscript tt (tr), the subscripts u 
and V indicate two different kinds of basis functions in 
each e symmetry orbitals and Si takes either -1-1 or — 1 
depending on the orientation of the oxygen octahedron 
surrounding the i-th V ion site and has different sign 
between the two V ions in each V ion pair. The ener- 



gies of these orbitals Sdii') can be expressed in terms of 
the octahedral Dq and the trigonal Dtrg crystal fields 
parameters: 

Ediai) - -4i?<7 + 2Arg/3, 




Fig. 2. Schematic drawing of the monoclinic distortion of V2O3 
in the AFI phase. The circles denote V ions and the ovals show 
the V ion pairs. The allows represent the directions of the tilting 
of the V ion pairs in the AFI phase. 



In the AFI phase, the crystal structure is distorted 
from the corundum structure to the monoclinic one. This 
distortion causes tilting of the V-V bond direction of ev- 
ery V ion pair from the z-axis toward the negative side of 
the y-axis and breaks three fold rotational symmetry (see 
Fig. 2). A low symmetry crystal field, which hybridizes 
the Cy and ai orbitals and also the Cy and e„ orbitals on 
the same site, is present in this phase. Although there 
are five independent parameters for this low symmetry 
crystal field C,ie^,ai), C^{el,eZ), a(e^,ai), C^{eZ,el) 
and Ci(eJ,e2), the most important are those related to 
the orbitals: Ci(e^,ai) and Ci(eJ^,eJ). Other three 
parameters involving orbitals are, for simplicity, all 
assumed to be zero. Note that there are relations be- 
tween these parameters for the two V ion sites in the pair: 
CA{eZ.ai) = CB{eZ,ai), CA{eZ,eZ) = -CB{eZ,eZ). 

From recent analysis of high-energy spectroscopy spec- 
tra, V2O3 is not simply ascribed to the Mott-Hubbard 
type insulator, where the 3d orbital is the main compo- 
nent of the valence-band top and the conduction-band 
bottom, but intermediate between the Mott-Hubbard 
and charge-transfer type insulators, where both the valence- 
band top and conduction-band bottom consist -of iiighly 
mixed states of the 3d and O 2p orbitals.ca'E3'E3 It 
should be emphasized that although the model we as- 
sume here does not contain the O 2p orbitals explicitly, 
the effects of the hybridization between the 3d and 2p 
orbitals can be partially included by using renormalized 
parameters. As will be discussed next section, the on- 
site effective Coulomb energy between the 3d electrons 
should be strongly renormalized by this d-p hybridiza- 
tion. 
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2. 2 Effect of d-p hybridization 

Before discussing the details of the ground state prop- 
erties obtained from the model, in this subsection, we 
will briefly consider the effects of hybridization between 
the V 3d and the nearest neighbor O 2p orbitals. To 
this end, here we conmare a single V^^ ion model with 
a VOe cluster model,0 in which the V 3d orbital and 
2p orbitals of the six oxygens surrounding the V ion oc- 
tahedrally are included. The parameters adopted in the 
cluster model are the effective on-site 3d-3d Coulomb en- 
ergy [/dd=4.0 eV, the charge transfer energy between the 
3d and 2p orbitals A = 4.0 eV, the Slater-Koster param- 
eters for the 3d-2p hybridization pda = —1.73 eV and 
pdn = 0.75 eV and the octahedral crystal field lODq — 
1.0 eV. These parameter values are deduced from the 
high-energy spectroscopy experiments.ES' The configura- 
tions 3d2+"L" (n = 0, 1, 2, 3) are assumed for the VOe^" 
cluster, where L denotes a hole in the O 2p orbital. For 
the Slater integrals of the 3d-3d multipolc interactions 
and F^^, the atomic Hartree-Fock values reduced to 
80% are used for the VOg cluster model. 

Table I. Eigen energies relative to the ground state energy for the 
lowest ten states calculated with the VOg~ cluster and the V'^"'" 
ion models are shown. The representation in the Oh symmetry 
and dominant electron occupation in the 3d? configuration of 
each eigen state (enclosed by the parentheses) are also listed. 



Eigen state Energy (eV) 
Cluster Ion 



3Ti 














1.40 


1.37 






1.42 


1.40 






1.70 


1.53 




(t2geg) 


2.68 


2.82 


lAi 


K) 


2.81 


2.84 


IT2 


it2geg) 


3.14 


3.10 


iTi 


(t2geg) 


3.39 


3.23 






3.57 


3.42 




(4) 


4.86 


4.77 



In Table I, the lowest ten eigen energies with respect 
to the ground state energy calculated with the VOe^^ 
cluster and the single V'^"'" ion models are listed. For 
simplicity, the 3d spin-orbit interaction is not considered 
in both the models. For the V^^ ion model, the Slater 
integrals F^^^ reduced to 90% of those for the cluster 
model and lODq — 1.7 eV are adopted for the V'^+ ion 
model. In the ground state of the VOg^" cluster model, 
the weights of the configurations 3d^, 3d^L, 3d'^L^ and 
3d5L3 are 41.8%, 44.6%, 12.5% and 1.1%, respectively. 
Although such strong 3d-2p hybridization, it is clearly 
seen in Table I that the energy level scheme of the VOg^^ 
cluster is well reproduced by the V'^"'" ion model with 
the renormalized parameters Fj^^ and lODq. Similar 
low laying level correspondence was found between the 
V06^°" cluster and {3d^) ion calculations and also 
between the VOg*" cluster and V^^ (3d^) ion calcula- 
tions. Because of this correspondence of the energy levels 
of the two, the V ion model would be good starting point 
to describe the electronic state of V2O3. 



The effective on-site Coulomb energy between "3d elec- 
trons" for the V ion model is greatly reduced compared 
to the model with the O 2p orbitals. The on-site ef- 
fective Coulomb energy between two different ''t2g or- 
bitals" U' can be estimated from the VOg cluster model 
by & = i;g(V0g«-)+£;g(V0gi°-)-2Sg(V0g^-), where 
Eg^s denote the ground state energies of the clusters with 
different charges. The on-site Coulomb energy between 
electrons on the same "i2g orbital" U and the on-site ex- 
change energy of the "t2g orbitals" J was obtained Ly 
using U ^ U' + 6F2 + 40Ei and J = 3F2 + 20F4H> 
The values obtained are U = 3.61 eV, U' = 2.31 eV 
and J — 0.65 eV. Similarly, the renormalized value for 
Udd was obtained by Udd ^U' + ff'a - 10^4=2.07 eV. 
These are the values which should be used in models in 
which the O 2p orbitals are not explicitly included. For 
comparison, these quantities for the bare t2g orbitals ob- 
tained with Udd = 4.0 eV and reduced to 80% from 
Hartree-Fock values are U = 5.60 eV, U' = 4.15 eV and 
J = 0.73 eV. 

Table II. Values of the renormalized parameters adopted in the 
V ion pair model. For the d-d hopping integrals and the crystal 
field parameters CA(e5},ai) and CA(eu,eJ), values for the AFI 
and PM phases are indicated. All values are units of eV. 









(AFI) 


(PM) 




2.07 


dda 


-0.735 


-0.793 


^dd 


7.29 


dd-K 


0.100 


0.109 


^dd 


4.57 


ddS 








lODq 


1.7 


CA(e:,ai) 


0.02 





Q 


0.023 


CA{eZ,eZ) 


-0.02 






Because of this strong d-p mixing, the hopping inte- 
grals dda, ddn and dd6 between the two V ion sites in 
the pair must be also renormalized. Effective d-d hopping 
matrix elements, where the influence of the d-p mixing is 
considered within the second-order perturbation theory, 
are 

- (^^- I ^/6W(pda)^ 

3^N^,\ A )' 

where 7 = u, w, A^^ = 1 + 4(pd7r)^/A and A^o- = 1 + 
6{pda)'^ / A. Using parameter values for pda, pdn and A 
deduced from the high-energy spectroscopy experiments 
mentioned before, it is found that values of the matri- 
ces other than V{ai,ai) is largely reduced from those 
without d-p hybridization. Such the effects of the d-p 
hybridization was included approximately by reducing 
the value of ddTr to 25% from the bare value. Harrison's 
values are adopted as for the bare values of the d-d hop- 
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ping integrals The renormalized parameters adopted 
in the V ion pair model are summarized in Table. II. 

2. 3 Properties of ground state 
In thi s subsection, we will discuss the ground sta te 
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Fig. 3. Trigonal crystal field -Dtrg dependence of the ai orbital 
occupation number per V ion in the ground state for the PM (a) 
and AFI (b) phases are shown. The vertical line is positioned 
at -Dtrg = 0.28 eV where the ground state gives the ratio e'^ai : 
= 1 : 2 deduced from the XAS experiment for the AFI 

phase. 



p rope rties of the V2O3, using the model described in 
§ 2T. In Fig. 3(a), the trigonal crystal field Z?trg depen- 
dence of the oi orbital occupation number for a V ion in 
the ground state for the PM phase is shown. Two dis- 



continuous decrease in the occupation number at -Dtrg = 
0.09 eV and -Dtrg = 0.34 eV are seen as -Dtrg increases. 
Because of the Hund's rule coupling, the spins of two 
electrons in each V ion site form S ^ 1 spin state in the 
ground state. This is very diffijrent result to the model 
proposed by Castellani et al.f3 where one electron-spin 
in the ai orbital couples to that of the other V ion of 
the pair forming a singlet spin state with the neighbor- 
ing site and 5=1/2 spin degree of freedom arising from 
the other electron in the the orbital remains in each V 
ion site. The difference in the spin state of the V ions is 
essentially caused by their choice of the on-site exchange 
interaction value J ~ 0.2 eV, which is less than one third 
of our value. For -Dtrg < 0.09 eV, where energy of the 
oi orbital is lower than that of the e'^ orbital, one of 
electrons is in the e'^ orbital and the other is in the ai 
orbital in both V ion sites (denoted by e^ai;e'^ai). On 
the other hand, for -Dtrg > 0.34 eV, two e'^ orbital are 
occupied in both V ion site (denoted by e^e^;e^e'^). In 
between these two region, the ground state is a superpo- 
sition of the e'^e'^;e'^ai and e'^ai;e'^e'^ states. 

For 0.09 < -Dtrg < 0.34, the ground state is doublet 
state and their wave functions \g±) in large U and lODq 



limit can be approximately described by 

l5+> = ^(|e;T«iT;<Te-T) + Kt^-t; <T«it))' 
Iff-) = ^{\el^a,^;el^el^) + \el ,el el ^a, ,)) , 

where t (i) denotes the magnetic quantum number of 
the electron spin s^ = 1/2 (s^ = —1/2) and |ej) — 
-i=(ip|e^) -i|ej)) with {el\L^\el) = =f1- These states 
nave the orbital and spin angular momenta along the z- 
axis and the values per V ion are {g± \Lz\g±) = and 
{g±\Sz |ff±) = ±1. Since this degeneracy of ground state 
is magnetic, there is no orbital degree of freedom. This 
is very different result to .the similar V ion pair model 
propose by Mila et al.!iBc3 where the orbital angular 
momentum is assumed to be small and thus the orbital 
degree of freedom still remains apart from the spin degree 
of freedom. The ground state wave functions of their 
model are written as 

\gr,Sj) = i=(|e-ai;e„e:) + |eX;e>i)) <S> = 2;Sj), 

where 7 specifies the two kinds of the e'^ orbitals u and v 
and IS""^ = 2; Sj) is the total -spin wavefunction, in which 
5* = 1 Hund's spins of the two V ions are coupled fer- 
romagnetically. The ground state wavefunctions of our 
model \g±) can be describe in terms of the wavefunctions 
I57; Sj) as 

l5±> = ■^{-i\gu;Sj = ±2)±\gv;Sj = ±2)). (2) 

The wavefunctions of Mila's ground state other than 
above two are lifted higher energies ranging up to ^ 
0.027 eV. 

It has been reported from the x-ray diffraction exper- 
iment that the distortion from the corundum structure 
to the monoclinic structure at the PM — > AFI transi- 
tion involves abrupt expansion of the V-V bond length 
of the V ion pair and rotation of the V-V bond direc- 
tion about 1.8° from the C/j-axis (within the y-z plane in 
Fig. 2).EP To take into account the effects of the lattice 
distortion, the values of dda and ddir are reduced from 
those for the PM phase using the V-V distanpa R depen- 
dence (ddm) oc proposed by HarrisorO^ and the 
low symmetry crystal field are considered as in Table II. 

In Fig. 3(b), the trigonal field -Dtrg dependence of the 
oi orbital occupation number in the ground state for the 
AFI phase is shown. The most different part compared 
to Fig. 3(a) is that the occupation number decreases 
continuously as -Dtrg increases around -Dtrg ~ 0.3 eV. 
This indicates that the boundary between the e'^e^; e'^e^ 
and e'^ ai, e'^ e'^ [e^e"; e^ai) states is obscure and around 
Arg ~ 0.3 eV, there is e'^e'^; e^e^ and e^ai; e'^e'^ {e^e""; e'^oi) 
mixed state. This can be understood as follows. Because 
of the symmetry lowering of the pair caused by the tilt- 
ing of the V-V bond direction from the trigonal axis, 
the ai and orbitals hybridized in each other. This 
leads to hybridization of the e^e'^\e^e^ and e'^ai\e^e'^ 
(e'^e^; e^ai) states around Dtrg 0.3 eV, where energies 
of these two states are almost degenerate. It is also seen 
in Fig. 3(b), the boundary is shifted to much lower -Dtrg 
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than that for the PM phase. This shift of the boundary 
is due to reduction in the V-V hopping integral values. 

If we assume -Dtrg ~ 0.3 eV, the abrupt change in the 
ai orbital occupation number is expected at the metal- 
insulator (M-I) transition. Such a large difference in the 
Oi orbital occupation number between the metal and in- 
sulating phases was indeed inferred from the measure- 
ment of thnJinear dichroism of the V 2p X-ray absorption 
spectrumJia^ Both the e'^ai to e'^e'^ ratio estimated from 
the measurement for the AFI phase e^ai : = 1:2 

and the PM phase e^ai : = 1:1 are well repro- 

duced with Z^trg = 0.28 eV. Since large change in the 
electronic state are caused by the lattice distortion, i.e., 
elongation and rotation of the V-V bond of the V ion 
pair, in this model, it is expected that the lattice distor- 
tion is strongly related to the metal-insulator transition. 
It is also important to mention that there is also large 
difference in the weight of the d^d}- and d^d"^ configu- 
rations between the ground states obtained for the two 
phases. While the weight is w(c?^rf^) — wid^S") — 3.5% 
for the AFI phase, the weight for the PM phase is 6.9%. 
Such large increase of the charge fluctuation from the 
AFI phase to the PM phase in this model again indicates 
strong connection between the monoclinic distortion and 
this met al-insulator transition 
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Fig. 4. Trigonal crystal field -Dtrg dependence of the (ey) or- 
bital occupation number n(eu) (n(e„)) per V ion in the ground 
state without (a) and with (b) the 3d spin-orbit interaction 
for the AFI phase are shown. The vertical line is positioned 
at -Dtrg = 0.28 eV, where the ground state gives the ratio 
e'^ai : e'^e'^ = 1 : 2 deduced from the XAS experiments for 
the AFI phase. 



To see the effects of the 3d spin-orbit interaction on 
the ground state wavefunction in the AFI phase much 
clearly, in Fig. 4, Dtrg dependence of the e„ and e„ orbital 
occupation numbers n(e„) and n(ev) per V ion without 
(a) and with (b) the 3d spin-orbit interaction are shown 
for the AFI phase. There is large difference between the 
two cases. For 0.09 < -Dtrg < 0.25, without the spin- 
orbit interaction, the occupation numbers are n{eu) ^ 1 



and n(ey) ~ 1/2 and this shows that the ground state 
is mainly dominated by the e^oi; e'^e'^ (e^ej; ej^ai) con- 
figuration. This also indicates that the monoclinic dis- 
tortion will induce the ferro V2 molecular orbital Afdfjr 
in the AFI phase, which is proposed by Mila et alS3\}3 
On the other hand, in the same parameter range, with 
the spin-orbit interaction, the occupation numbers are 
n{eu) ~ n{ev) ^3/4 and both two kinds of the config- 
urations ejai;ejej and ejai;ejej are almost in equal 
weight. As was discussed before, this is because the 
strong coupling between the spin and orbital degener- 
acy and there are no orbital degree of freedom which is 
independent to the spin degree of freedom in this ground 
state and thus we do not expect usual Jahn- Teller insta- 
bility. Around -Dtrg ~ 0.3 eV, where the e^e^; e^e^ state 
is mixed with the e^ai\e^e^ [e^ \ ai) state, the sit- 
uation is similar: difference between n(eu) and n{ey) is 
smaller with the spin-orbit interaction than that without 
the interaction. 

The ground state property with the spin-orbit inter- 
action above can be understood more clearly with the 
aide of an effective Hamiltonian within the limit of large 
U and lO-Dq. In this limit, the ground state is well de- 
scribed within restricted basis functions \ge) = |e^e^;ej^e^) 
and I57) = ■^{\e^^ai;elel) + \eleZ-e-^^ai)) (7 = u,v). 
For convenience, here, L' = 1 pseudo orbital angular 
momentum operator L' is introduced, where the eigen 
functions L'^\g m) = m\g m) of L'^ are defined hy \g m = 
0) = -\ge) and \g m = ±1) = {-i\gu) ± \gv))/y/2. The 
effective Hamiltonian can be expressed using L' and the 
S"^ — 2 total spin operator 5'"'^ of the V ion pair as 



+ V2CA{e:,a,)iL',L'y + L'yL',] 



(3) 



where A = C3d/4and A = Dt,f,~{idday + {dd7r/3f}/{U'- 
J). Note that there are relations between the real total 
orbital angular momentum operator L and L' within 
this restricted basis: — —V2L'^, Ly ~ —\f2L'y and 
Lz — —L'^. If CA(eJ,ai) = 0, the system has uniaxial 
anisotropy and when A < — 2A, the ,z-axis is the mag- 
netic easy axis and the ground state is \g±) = |gm = 
ibl)!^^ = ±2), which is identical to that in eq. (||). 
On the other hand, for A > — 2A, the x-y plane is the 
magnetic easy plane. With finite CA(eJ^,ai), the term 
[L'^L'y + L'yL'^) causes hybridization between \ge) and 
\gu) and this rotate the magnetic easy axis from the 
z-axis toward the y-axis. The ground state can be re- 
garded as magnetic doublet state having their magnetic 
moments along this easy axis with opposite directions in 
each other, although breaking of the rotational symme- 
try around this easy axis causes a small energy splitting 
< 0.5 meV between these lowest two eigen states. 

2.^ Magnetic anisotropy 

As was discussed in the previous subsection the admix- 
ture of the e"e"\e^e^ state in the ground state lead to 
rotation of the magnetic easy axis. To study such mag- 
netic properties of the ground state, a molecular field Hm 
with a canting angle 9m with respect to the z-axis in the 
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Fig. 5. Sizes of the spin Ms, orbital A/l and total A/t magnetic 
moments (a) and the canting angle of the total magnetic moment 
from the z-axis 8m (b) as a function of the trigonal crystal field 
Dtrg- 



y-z plane was further applied in the V ion pair model: 
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In the calculation, the molecular field with the magni- 
tude Hm = 0.01 eV was adopted and 9m was chosen 
so as to give the lowest ground state energy. The mag- 
netic moment is, therefore, parallel to both the magnetic 
easy axis and the molecular field. Note that the values 
and direction of the spin and orbital magnetic moments 
obtained here are not sensitive to the magnitude of the 
molecular field. 

In Fig. 5, the sizes of the spin A/g, orbital A/l and 
total Afx magnetic moments per V ion (a) and 9m (b) 
are shown as a function of i^trg- Both the directions of 
the magnetic moment and the V ion pair are in the y-z 
plane. For the e^ai^e^e^ {e^e^;e'^ai) ground state, the 
direction of the magnetic moment is parallel to the z-axis 
(the corundum c^-axis). However, as the weight of the 
gTTgTT. gTTgTT configuratiou in the ground state increases, 
the magnetic moment is canted toward the y-axis and fi- 
nally reaches to the y-axis with a pure e^e'^; e'^e'^ ground 
state. Note that the V-V bond of the V ion pair is tilted 
toward opposite direction to the magnetic moment, i.e., 
toward the negative direction of the y-axis. The large 
orbital magnetic moment A/l ~ 0.7/iB is present in the 
e'^ai;e^e^ {e^e^;e^ai) ground state. The orbital mag- 
netic moment is about 40% larger than that expected 
from the ground state wavefunction (see eq. (2) ) in the 
limit of large U and lODq. This deviation is mainly 
due to the presence of the small weight of configurations 
with occupied e"' orbitals arising from the 3d-3d multi- 
pole interaction. While Af^ remains almost unchanged 
Ms ~ 2/iB as the weight of the e^e^;e^e^ configura- 
tion increases. Ml decreases rapidly. At Dtrg — 0.28 eV, 



Mt ^ 1-2/iB and 9m = 51° and these well agree with the 
values obtained in the neutron diffraction experiments 
Mt ^ 1.2/iB and 9m ~ 71°B Strictly speaking, the 
magnetic moments of the two V ions are slightly tipped 
out from the y-z plane toward opposite direction in each 
other. Although the x-component of the magnetic mo- 
ments is an order of ^ 0.01/iB, the fact is important for 
the interpretation of the RXS measurements as will be 
discussed in 



5.3 



§3. Electron-Lattice Interaction 

As was discussed in § the pure e^oi; e'^e^ (e^e^; e^ Qi) 




Fig. 6. In (a), tilting angle from the cji-axis as a function of the 
V-V distance of the V ion pair R is shown for G = 7.85, 15.7 and 
31.4 meV/ (°). -Rafi pi denotes the experimentally observed V-V 
distance of the V ion pair in both the AFI and PI phases whereas 
RpM shows that in the PM phase. In (b), the ai, Su and e„ 
orbital occupation numbers calculated with G = 15.7 meV/(°) 
are shown. 



ground state has only spin degree of freedom and no or- 
bital degree of freedom is left behind. One would not, 
therefore, expect usual Jahn- Teller lattice instability in 
this system. However, in the AFI phase, V2O3 exhibits 
the monoclinic lattice distortion and this indicates that 
there must be some other kind of lattice instability in 
this system. Here, it will be shown that such instabil- 
ity indeed appears when the energy of the e'^ai;e'^e^ 
{e^e^;e^ai) and ; states are almost degener- 
ate. 

In order to discuss interaction between the lattice de- 
formation and the electronic state, here, we consider one 
V ion pair in the crystal and in addition to the energy 



of the electronic state discussed in § 2.1, the local lattice 



8 



Author Name 



deformation energy caused by tilting of the V ion pair is 
included. The energy of the V ion pair as a function of 
the polar 6 and azimuthal angles of the tilting direction 
would be written in the form: 



Etot{e, if) = E,x,,{e, p) + Y^'' + FO-" sm3ip + ... ,(4) 

where the first term denotes the ground state energy of 

and the second and 



2.1 



the Hamiltonian described in 
third terms are the local lattice distortion energies. 

Instead of including F0^ sin 3(^9 and much higher order 
trigonal potential terms of the lattice deformation, for 
simplicity, we consider the lowest term only and as- 
sume that tilting occurs only to the restrict directions of 

= — -I, — -fi^, where displacements of each V ion 
of the pair is parallel to one of the three V-V bond in 
the honeycomb layer. As a result, the V-V bond along 
this direction is elongated and those in the other two 
directions contract. Note that this tilting direction is 
consistent with that observed in the monoclinic lattice 
distortion in the AFI phase. In such lattice distortion, 
both the ej^ and e'^ orbitals can be hybridized with the ai 
orbital and corresponding low symmetry crystal field pa- 
rameters, which are linear functions in 6 for small lattice 
distortion, can be written in the form: 

CAieZ,ai) ^ ^cosif, CA(eJ^, ai) = ^ sin (5) 

In the calculation. Dug = 2.8 eV was adopted and the 
V-V distance R dependence for the d-d hopping integrals 
dda, ddn ex. was assumed. 

In Fig. 6, the tilting angle at the minimum of i?tot as a 
function of the V-V distance of the V ion pair R is shown 
for three different crystal fields G = 7.85 (open circles), 
15.7 (closed circles) and 31.4 meV/(°) (triangles). The 
coefficient K was chosen so that maximum tilting angle 
is 1.8° and those for G = 7.85, 15.7 and 31.4 meV/(°) 
are K = 2.60, 6.50 and 15.7 meV/(°)2, respectively. The 
V-V distances i?AFi,Pi and RpM are the experimentally 
observed values in the insulating phases and in the PM 
phase, respectively .B The rotation of the V ion pair oc- 
curs only when the energy of the two kinds of the states 
e'^oi; e'^e^ {e^e^; e^ai) and e'^e'^; e^e" are nearly degen- 

o 

erate and takes the maximum value around R = 2.750 A, 
where the two kinds of the states are just degenerate. 
Note that the energy difference between the two can be 
given by - Arg - {[ddaf -f [ddn /if} / {U' - J). 

This instability toward the lattice distortion can be 
explained by the similar mechanism to the Jahn- Teller 
effects. If we include the lattice distortion energy, the 
effective Hamiltonian eq. (3) generalized for the V ion 
pair with an arbitrary tilting angle {6, ip) can be written 
as 



,ff = ~X{SjL', + V2{SjL'^ + 5Jl;)} + AL'l 

+ Ge{cos p{L',L', + L',L',) -f sin p{L',L'y + L'^L',)} 
K 



-9^ +F0^ sinSip 



(6) 



When the energies of the these two kind of states \ge) 
and \gj) (7 = u,v) are degenerate or nearly degenerate. 



the term linear in 9, which couples the states \ge) and 
\g"f), causes rotation of the V ion pair in the same way 
as the Jahn- Teller effects. The condition for instability 
toward tilting is approximately obtained for small A and 
\F\ as |A| < 2G'^/K. 

Since the energy gain from this electron-lattice inter- 
action is maximum at the V-V distance where the ener- 
gies of the two kinds of the states are nearly degenerate, 
if this mechanism is the cause of the monoclinic distor- 
tion in the AFI phase, the V-V distance of the pair in 
this phase should be at this distance. Indeed, the value 

o 

R — 2.755 A giving the maximum tilting angle is very 

o 

near to the value i?AFi,Pi = 2.745 A, where the V 2p 
XAS experimental results on the e^e'^ and e'^ai config- 
uration ratip_for the AFI phase e^ai : e^e^ = 1 : 2 is 
reproduced. EfP The monoclinic lattice distortion in the 
AFI phase thus can be understood by this mechanism. 
On the other hand, the V-V distance is much shorter in 
the PM phase i?PM - 2.700 A and if |G| < 16 meV, the 
lattice instability does not occur at this distance. In this 
situation, pure e^ai;e^e" [e^e^^e^ai) state is expected 
in this phase. This also agrees with the V 2p XAS ex- 
periment, where e^ai : e^e^ = 1 : 1 is obtained in the 
PM phase. 

Although the experimentally observed V-V distance is 
the same to that for the AFI phase, in the PI phase, no 
monoclinic lattice distortion has been found. This can be 
explain as follows. In the PI phase, each V ion pair is also 
tilted from the c/j-axis similar to the AFI phase. How- 
ever, their tilting directions are disordered and fluctuate 
among the three stable directions (^ = — ^, —^zt^, 
resulting in no monoclinic lattice distortion. This as- 
sumption is consistent with the V 2p XAS experimental 
results, where the weight of the e'^e^ configuration is 
larger than that of the e^ai configuration both in the PI 
and AFI phases in their initial state, indicating admix- 
ture of the e'^ e'^ ; e'^ e'^ configuration. If this assumption 
is correct, the AFI PI transition is accompanied by 
order to disorder transition in the tilting directions of the 

V ion pairs. As will be discussed in the next section, this 
ordering of the tilting direction in the AFI phase can be 
explained by a model in which this local lattice distor- 
tion and the magnetic exchange interaction between the 

V ion pairs are considered. 

From above discussion, it is expected that there are 
at least two local minima for the free energy of the real 
system as a function of the V-V distance of the pair. 

o 

One is positioned at RpM ^ 2.700 A, where the energy 
difference of the two kinds of the state e'^ ai; e'^ e'^ and 
gTTgTT. gTTgTT jg Jargc aud thus no electron-lattice coupling 
takes place, and the other one is located at i?AFi,Pi ~ 

o 

2.745 A, where energy of the two kinds of the states 
is degenerate and Jahn- Teller like lattice distortion is 
the cause of the energy lowering at this V-V distance. 
The expansion of the V-V distance in the AFI and PI 
phases as compared to that in the PM phase can be 
understood as a consequence of the large change in the 
3c? orbital occupation owing to the tilting of the V ion 
pair. Since mixture of the e'^e'^; e'^e'^ state in the ground 
state reduces electron occupation of the ai orbital, which 
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point toward the direction of the V-V bond, the nuclear 
charges of the V ions of the pair is much less screened 
by the 3d electron in the V-V bond direction, resulting 
in elongation of the V-V bond. 

Since both in the AFI and PI phases, the electron- 
lattice coupling strongly influences the orbital occupa- 
tion and charge fluctuation, the PM AFI and PM 
PI transitions can not be regarded as usual Mott tran- 
sition. In S 2.3, it is found that when the lattice distor- 



tion occurs, charge fluctuation between the two V ions 
is strongly suppressed. The expansion of the V-V dis- 
tance and resultant reduction of the d-d hybridization 
strength due to this electron-lattice coupling in the insu- 
lating phases is probably the main cause of this metal- 
insulator transition. 

§4. Magnetic Interaction between V Ion Pairs 

In the previous section, we have discussed about the 




T; =C0 
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Fig. 7. Mapping of V ion pairs in the corundum lattice on the 
simple cubic lattice (a) and three possible directions correspond- 
ing to Ti = Lui, i02 and uj^ of the tilting of the V ion pair from the 
Ch-axis (b) are drawn. On left hand side of Fig. (a), the ovals 
show the V ion pairs and the allows represent the spin direc- 
tions of the V ions and on right hand side, the directions of the 
pseudo spin (Ti = ±1 for V ion pairs are indicated by the arrows. 
In Fig. (b), the circles and arrows represent the V ions and their 
displacement due to the tilting of the V ion pair, respectively. 



local lattice distortion effects on the electronic states 
within individual V ion pairs. Here, we will consider 
about the exchange interaction between the V ion pairs 
in its relation with the local lattice distortion in order to 
discuss the AFI PI transition observed in Vi-a^CrajOa. 
The transition is accompanied by the monoclinic to corun- 
dum lattice structure transition and in the ultrasonic 
wave velocities measurements, the C44 elastic constant 
was estimated to soften almost at the transition tem- 
perature, shxming the nearly second order nature of the 
transitionJlJ' t3 Contrary to this, in the AFI to PM 
transition such softening of C44 is not observed and the 



transition is of the first order. In the magnetic thermal- 
neutron scattering experiments, abrupt change in short 
range spin correlaticm_£unction at the transition temper- 
ature was observed.oll3) 

In the AFI phase, the monoclinic lattice distortion and 
the magnetism are connected, forming unique AF mag- 
netic order. In the monoclinic distortion, all the V ion 
pairs, which is parallel to the c/j-axis in the coxundum 
structure are tilted ^ 1.8° in the same direction.B With 
this displacement of the V ions, the honeycomb lattice, 
which has three hold rotational symmetry in the corun- 
dum structure, is distorted and V-V bonds in one of the 
three directions in the layer are elongated and those in 
the other two directions contract (see Fig. 2). The mag- 
netic moments between two V ions with the elongated 
V-V bond on the honeycomb lattice are ferromagneti- 
cally coupled and those with the contracting V-V bond 
are antifcrromagnetically coupled. 

In the thmdes proposed by Castellani et alx3 and 
Mila et al.l^3i3 the orbital ordering is an essential in- 
gredient for the formation of this unique AF magnetic 
order in V2O3, where exchange coupling constant be- 
tween the spins of V ions depends on the orbital oc- 
cupation of the V ions. Rice has pointed out that the 
abrupt change in the spin correlations at the AFI PI 
transition and the short-range spin fluctuations in the 
PI phase can be explained in the presence of the orbital 
degree of freedom. Contrary to this, Ezhov et a/O 
performed LDA+J7 band structure calculations and ob- 
tained a solution with the unique AF magnetic order 
as the lowest total energy magnetic state, assuming the 
monoclinic lattice distortion. Since the orbital occupa- 
tion in each V ion is e'^e'^, there are no orbital order 
in this AF magnetic state. In the present theory, since 
strong coupling between the spin and orbital degrees of 
freedom due to the spin-orbit interaction, orbital order 
does not exist in the ground state. Here, it is assumed 
that the coupling constant of the exchange interaction is 
a function of the V-V bond length and not depends on 
the orbital occupation of the V ions. 

4-1 Effective spin-lattice Hamiltonian 

The effective spin Hamiltonian considered here con- 
sists of the exchange interaction between the nearest 
neighbors of the V ion pairs and the local lattice distor- 
tion energy in each V ion pair. Each V ion pair has six 
nearest neighbors and as depicted in Fig. 7, the network 
of the V ion pairs can be mapped onto the simple cubic 
lattice. Two spin directions of the doublet ground state 
on the i-th V ion pair is projected onto Ising like pseudo- 
spin with (7i = ±1 at the corresponding lattice point i of 
this simple cubic lattice. The coupling constant of the 
exchange interaction is assumed to be a function of the 
displacement of the V-V bond length due to tilting of 
the V ion pairs. The effective spin Hamiltonian then can 
be written as 

ihj) * 
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where the first term represents the exchange interaction 

between the pseudo-spins of the V ion pairs and the sec- 
ond term describes the local lattice distortion energy in 
each V ion pairs. 0i >0 represents tilting angle of the i- 
th V ion pair and two dimensional unit vector Tj denotes 
the direction of the tilting projected on the honeycomb 
layer. Here, we assume that each V ion pair is tilted in 
such a direction that each V ion in both end is parallel to 
one of the three V-V bond in the honeycomb layer and 
elongates this V-V bond. Tj is, therefore, in one of these 
three directions: u^i = (^/3/2,l/2), LO2 = (-%/3/2, 1/2) 
or = (0,-1). lij — LJi, (jJ2 or ui^ indicates the direc- 
tion of the V-V bond between the i- and j-th V ion pairs. 
Thus, the term (0jTj + OjTj) ■ lij is proportional to the 
elongation of the V-V distance between the i- and j-th V 
ion pairs. 

The AF magnetic order observed in V2O3 corresponds 
to the C-typo AF magnetic order in this model, where 
the nearest neighbor spins along one of the lattice line 
in the simple cubic lattice are coupled ferromagnetically 
and those along the other two are antiferromagnetically 
coupled (see Fig. 7 (a)). For this AF magnetic order, the 
energy per V ion pair is easily obtained as a function of 
a uniform distortion 9 = 9i 

E{e) = -\J\C - (G + A\J\)e + 



and the energy has minimum at ^mii 

1 (G- 



G+4\J\ 
K 



E{9„ 



-\J\C 



2 K 

For this AF magnetic state being stable, the exchange 
coupling constant J(2^i„in — C) for the long V-V length 
should be positive and that for the short V-V length 
— J{0min + C) should bc negative. Thus at least the con- 
ditions J < and 2^ini„ > C > — ^min must be satisfied 
for this AF magnetic state being the ground state. 

4.2 AFI to PI phase transition 

To study finite temperature properties of this effective 
Hamiltonian, Monte Carlo simulations were performed 
with a 24 X 24 X 24 periodic lattice for various values 
of G. Figure 8 shows the phase diagram as a function 
of G and kBT/\J\ obtained from the calculation. The 
adopted value of G is G = 15.7 meV/(°) and from the ex- 
perimental facts 9 

min 

1.8° and Tn ~ 155 K, the other 
parameter values are estimated to be J ^ —4.0 meV 
/(°), K/\J\ - 3.9 and G/\J\ - 3.0. The solid 

and dashed lines represent the phase boundaries calcu- 
lated with {G/\J\ = 3) and without (G/|J| = 0) the 
electron-lattice coupling, respectively. In the range of 
—0.3° < G < 2.25° as shown in the figure, the low tem- 
perature phase is that with the C-type AF magnetic or- 
der and the V ion pairs arc tilted and ordered in the same 
direction in this phase (denoted by CO). This phase cor- 
responds to the AFI phase in the real material. The 
high temperature phase is paramagnetic and the tilting 
direction of the V ion pairs are disordered (denoted by 
PD) and this phase corresponds to the PI phase in the 
real material. For C > 1.1°, there is an intermediate 
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Fig. 8. Phase diagram of the effective spin Ifainiltonian for 
GI\J\ = 3 (soUd line) and G/\J\ = (dashed Une). CO de- 
notes the C-type AF magnetic order with ordered tilting direc- 
tion of V ion pair. GD and PD are the G-type AF magnetic 
order and paramagnetic order with disordered tilting direction, 
respectively. 



temperature phase with the G-typc AF magnetic order, 
where all the nearest neighbor spins are antiferromagnet- 
ically coupled, with disordered tilting directions of the V 
ion pairs (denoted by GD). The transition from the CO 
phase to the PD phase is of the second order whereas the 
CO GD and GD PD phase transitions are of the 
first order. It is also seen in the figure that the CO phase 
is further stabilized with the electron-lattice interaction 
(solid line) than without it (dashed line), showing impor- 
tance of the interaction for the formation of the C-type 
AF magnetic order. 
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Fig. 9. Order parameter for the tilting of the V ion pairs (SiTj ■ 
Word) the tilting angle as a function of temperature 
calculated with C = 0.9° (solid hne) and C = 1.8° (dashed 
line). 



In Fig. 9, the order parameter of the tilting of V ion 
pairs {9iTi ■ Word) and the tilting angle {9i) as a function 
of temperature are shown for G = 0.9° (solid line) and 
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KTI\J\ (deg.) 

Fig. 10. I/xStCT) is plotted with error bars for C = 0.9° as a 
function of temperature. 



C = 1.8° (dashed line), where Word is one of the three 
possible tilting directions ijJi {I = 1,2,3) in which the V 
ion pairs are ordered. It is seen in the figure that the CO 
— > PD transition is of the second order and the CO GD 
transition is of the first order. The tilting angle {9i) in 
the PD phase is much larger than that in the GD phase, 
showing stronger dynamic lattice distortion in the PD 
phase. The order parameter of the tilting {9iTi ■ ujord) is 
expected to be approximately proportional to the degree 
of the monoclinic lattice distortion corresponding to the 
strain 64 (see Fig. 2). The reciprocal of the quantity de- 
fined by xeAT) = mT,-0Jori? ~ {{0,t, ■ uj,,df))lk^T 
is therefore proportional to the elastic stiffness constant 
C44. In Fig. 10(c), l/xeT{T) is plotted as a function of 
temperature. Both above and below Tn, I/xst decreases 
as temperature approaches to Tn- Since the transition is 
of the second order, the value of I/xst is zero just at Tn. 
This temperature dependence of 1/ xbt is consistent with 
the behavior QLC'44 observed in ultrasonic wave velocity 
measupancntstlP and also in neutron scattering exper- 
imentsO on Cr doped V2O3, where C44 softens with 
decreasing temperature in the PI phase and approaches 
zero near the PI — > AFI transition temperature. 

In Table III, spatial correlation functions (<JiCrj) and 
{OiTi ■ OjTj) between the first, second and third near- 
est neighbors along a lattice line are listed. Those for 
temperatures above and below the CO PD transi- 
tion calculated with C = 0.9° are shown and for com- 
parison, those for the CO GD transition calculated 
with C = 1.8° are also listed. In the CO phase, those 
along two different lattice lines, parallel and perpendic- 
ular to the ferromagnetic plane of the C-type AF order, 
are shown. In the PD phase, there is short range G-type 
AF spin correlation, which is consisterttjAfith recent neu- 
tron scattering experimental results. la'M In addition to 
this, short range order of tilting directions is also present 
in this phase. This sudden change in the spin correla- 
tion, which is unusual for order-disorder magnetic transi- 
tion can be explained as follows. In the PD phase, spins 



Table III. Spatial correlation functions {oitTj) and {diTi ■ 8jTj) 
between the first, second and third nearest neighbors along a lat- 
tice line are shown above and below the CO — > PD and CO — > 
GD phase transitions. In the CO phase, those along two differ- 
ent lattice lines, parallel and perpendicular to the ferromagnetic 
plane of the C-type AF order, are listed. 

C = 0.9°, K/\J\ = 3.9 {°)-i, G/\J\ = 3 
kBT/\J\ = 3.65° (CO phase) 

(cTiO-j) {OiTi ■ 9jTj) 

1st 0.678 -0.735 1.21 1.16 
2nd 0.670 0.702 1.07 1.09 
3rd 0.681 -0.691 1.03 1.05 

kBT/\J\ = 3.9° (PD phase) 

{(TiO-j) (diTi ■ 9jTj) 

1st -0.191 0.170 

2nd 0.141 0.032 

3rd -0.042 0.005 

C = 1.8°, K/\J\ = 3.9 (°)-i, G/\J\ = 3 
fcaT/IJI = 3.3° (CO phase) 

1st 0.964 -0.969 2.02 2.12 
2nd 0.965 0.966 2.10 2.04 
3rd 0.964 -0.965 2.07 2.04 

fegT/l J| = 3.8° (GD phase) 

(cTiO-j) {OiTi ■ 9jTj) 

1st -0.957 -0.001 

2nd 0.955 0.007 

3rd -0.955 -0.002 



and the local lattice deformation are still dynamically 
coupled and short range C-type AF order with local lat- 
tice distortion survives. However, since the short range 
C-type AF order with three different directions of the 
ferromagnetic plane are coexistence, we find, in average, 
the G-type AF spin correlation as seen in Table III. In 
this way, the abrupt change in the spin correlation func- 
tions observed at the AFI — > PI transition in the neutron 
scattering experiments is explained as a result of the si- 
multaneous order to disorder transition of both the tilt- 
ing directions and spins of the V ion pairs. On the other 
hand, in the GD phase, since lattice distortion is not re- 
quired for the stabilization of the G-type AF spin order, 
the correlation functions between the tilting directions 
{9iTi ■ OjTj) is almost zero. This is in contrast to the 
G-type AF spin correlation in the PD phase originating 
from the short range C-type AF order with the dynamic 
spin-lattice coupling. 

§5. Resonant X-ray Scattering Analysis 

Resonant x-ray scattering (RXS) has received much 
attention recently, as apteciiiiiaue for extracting direct in- 
formation on magneticHEjEa' and orbitalE^ order. Pos- 
sibility of observing the orbital ordering ii|i_V203 by RXS 
have been pointed out by Fabrizio et a/.EZP They made 
predictions on the light polarization and the azimuthal 
angle of the scattering plane dependence of the resonant 
Bragg reflection intensities, ^assuming the orbital order 
proposecLby Castellani et alS3 The RXS measurements 
on V2O3CJ) and its Cr doped alloysEil' have been per- 
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formed by Paolasini et al. They insisted that in the AFI 
phase, there are resonant Bragg reflection peaks arising 
from the orbital ordering in addition to those caused by 
the AF magnetic ordering. On the other hand, Lovesey 
and Knight claimed that all these Bragg peaks can be 
attributed to the AF magnetic ordering and are associ- 
ated with the orbitaLpagnetic moment and the octupole 
moment of a V ion.E3 In this section, discussion will be 
mainly focused on these resonant Bragg peaks, especially 
(111) reflection peak, whose origin is still controversial. 
It will be shown that the light polarization and the az- 
imuthal angle dependence of the (111) Bragg reflection 
intensity can be explain within the V ion pair model and 
is attributed to not the orbital ordering but the antiferro 
(AF) magnetic ordering. 

5.1 Is — > Ap dipole X-ray scattering process 




Fig. 11. Antiferromagnetic and crystal structures of V2O3 in the 
AFI phase. The white and black circles denote V ions and the 
arrows indicate direction of their magnetic moments. In (a), 
the monoclinic unit cell and the magnetic structure are shown. 
The magnetic moments are ferrromagnetically coupled in each 
monoclinic (010) layer, altering their direction between adjacent 
layers. The V ion pair and neighboring six V ions are depicted in 
(b). In (c) and (d), V ions belonging to two different (010) layers 
are separately depicted with surrounding oxygen octahcdra. 



^-th V ion site. 

Here, the 3x3 tensor F^uj) is defined by 

F^{io)^q^{g\P ^ 



Eg+Lu-Hr,, + ir,„/2 



(8) 



where \g) and Eg arc the ground state wave functions 
and its energy, and Fm are the Hamiltonian and a 
life-time term for the intermediate state reached by the 
Is — > 4p transition and i-' = er is the electric dipole 
moment operator. Since delocalized feature of the Ap 
orbital, the intermediate-state-multiplet dependence of 
RXS is not expected to be important. Thus, here we will 
consider the integrated value of F^{uj) over the photon 
energy lo: 



F' = 



FlAu:)duj. 



(9) 



Note that F^,,, F^^^ 



and FI 2 correspond to the Is ^ Ap 
dipole transition probabilities integrated by lo at site I 
with X-, y- and z- linearly polarized light, respectively. 

Before discussing rather complicate RXS in V2O3, first 
of all, let us consider the most simplest antiferromagnet, 
where the unit cell with one magnetic ion site in the ab- 
sence of magnetic order is doubled by the presence of the 
antiferromagnetic order and two magnetic sublattices A 
and B with opposite direction of the magnetic moment 
in each other are assumed. In this case, for an ion 
site belongs to the sublattice A (here we denote it by 
F^) can be obtained from that belongs to the sublattice 
B (denoted by F^) by replacing the ground state wave- 
function \g) in eq. (||) with 6T|g), where represents 
the time reversal operator and T denotes the transla- 
tion operator, which transforms ion sites belong to the 
sublattice A into those belong to the sublattice B. This 
immediately leads to the relation: 

F^ = F^* . 

For resonant Bragg reflection corresponding to the anti- 
ferromagnetic order, the integrated scattering amplitude 
/dd is 

/dd = Y Y^iKj - 

= iiV{(e.4 - ej,4)ImF4 + {sye', - e,e')lmF^^^ 



+ (e.e^ - e^e',)ImF^J, 



(10) 



where N denotes the number of the magnetic ion sites. 



ImF^y can be written as 



In this subsection, we will discuss the Is Ap dipole 
process, which is dominant process in the RXS at the V 
i^-edge absorption region. The x-ray scattering ampli- 
tude /dd i^) is expressed as 

— m-q'yR.^^'* .F\co)e), 



/dd(w) 



(7) 



where q and e are the wave vector and the polarization 
vector for the the incoming photon, q' and e' are those 
for the outgoing photon and Ri denotes position of the 



ImF^ 



2i 



t(F^i-FAi), 



where F^i are integrated transition probabilities of the 
Is — *■ Ap dipole transition with plus and minus helicity 
light injected along the z-axis for an ion site belongs to 
the sublattice A. Carra et al. have proposed a sum rule 
for the core-level x-rajj-absorption spectra with the cir- 
cularly polarized light .123^ From the sum rule, it is easy 
to find that the quantity F^-^ — F\ is proportional to 
the z-component of the orbital moment of the Ap orbital 
(L^P) for a magnetic ion site belongs to the sublattice 
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A in the ground state. Similarly, the quantities ImF.^^ 
and Im.F^^ are proportional to the x- and y-components 
of the 4p orbital angular momenta {L*p) and (Ly^), re- 
spectively. The integrated amplitude /dd is, therefore, 
proportional to (e x e') • {L'^^). 

In Fig. 11, the magnetic and crystal structures in the 
AFI phase are shown. The black and white circles are 
the V ion sites and the arrows indicate the direction of 
their magnetic moments. There are eight different V 
ion sites in the monoclinic unit cell shown in Fig. 11(a). 
In (b) and (c), two different (010) magnetic layers are 
separately drawn, and the eight V ion sites are labeled 
A to H. The integrated scattering amplitude /dd(/i, k, I) 
of RXS at a {h, k, I) Bragg reflection can be written in 
the form: 



fT,T>{h,k,l) 



+ {F^j + (-)"+'=+'i^5}e- 



(11) 



where a 47r(/ix + lz). x -0.3438 ^nd z ~0.2991 are 

the position parameters of the V ions.cP If we assume the 
antiferromagnetic structure above, the system is invari- 
ant under the symmetry operation {6| ^ (ctm + bm + Cm)} 
and by this operation, the V ion sites A, B, C and D are 
transformed into the sites F, E, H and G, respectively 
and as we have discussed above, relations = F^* , 
pB ^ pE*^ pC ^ pu* pD ^ pG* satisfied. 

Similarly, the V ion sites A, C, E and G are transformed 
into the sites B, D, F and H by the 12/ a space group 
symmetry operation {cThl^im} and relations F^ — F^, 
pC ^ pD pE ^ pF j^j^^ ^ j^j.g obtained. Using 

above relations among _F''s and also the fact that the 
operation Ch transforms {Sx^Sy, £z) into {£x,£y, —£z), we 
obtain the scattering amplitude for the reflections with 
h =even and h + k + 1 =odd as 

fuD{h,k,l) = 

^ {(£,£; - e.e'y)ilmFy^^, + e'^ImF^) 

+ (e,4 - e,4)(Imi^,^, + e^ImF^J}. (12) 

On the other hand, the scattering amplitude for the re- 
flections with h =odd and h + k + l =odd can be written 
as 



fBi)ih,kJ) 
2 \^v^z 



.4)(ImF, 



A 



nmF^J. (13) 



It is obvious from eqs. (12) and (13) that the ampli- 
tudes are zero for the a-a scattering, where the polariza- 
tion vectors e and e' are perpendicular to the scattering 
plane and thus are parallel to each other. This is consis- 
tent with the experiments, where only vanishingly small 
intensities of both the (111) aad (221) reflections for the 
a-a scattering are observed.© On the other hand, the 
(T-TT scattering, where one of the polarization vectors is 
perpendicular to the scattering plane and the other is 



in the plane, is not forbidden scattering and indeed, for 
the magnetic (221) Bragg reflection, strong RXS are ob- 
served in the experiment. The experimental ct-tt scat- 
tering at the (111) Bragg reflection, however, has only 
negligible small intensity in the Is — > Ap excitation en- 
ergy region. This fact can be explain as follows. Since 
the direction of the magnetic moment observed in the 
neutron experiment is almost in the a„i-bm plane, the 
x-component of the Ap orbital angular momentum (L^^) , 
which is perpendicular to the a,n-h„i plane (see Fig. 11), 
is expected to be zero or small. Because eq. (13) only 
contains Fy^ and F^.,, which is proportional to (i^^), 
the intensity should be almost zero. 

5.2 Interference of Is — s- 4p and Is — s- 3d processes 

We have discussed the Is ^ 4p dipole RXS. Only 
the dipole process, the scattering amplitude at the (111) 
Bragg reflection should be almost zero and thus the Is ^ 
3d quadrupole transition expected to be important to 
explain the sharp peak observed at the prethreshold re- 
gion. Since V^^ ion sites do not have inversion sym- 
metry in V2O3, the Ap and 3d orbitals on the same 
site can be hybridized in each other and this causes in- 
terference between the Is — > 4p dipole trajasition and 
Is — !■ 3d quadrupole transition processes.EZP Because 
pure Is 3d scattering amplitude is expected to be 
much smaller than the interference scattering amplitude, 
here, only this interference scattering process is consid- 
ered. 

The crystal field which causes the hybridization be- 
tween the Ap and 3d orbitals in the same V ion site would 
be written as 

V'{r) = ± Ai,orYi,o(r!) ± As^or^Ys^oin) 

±A3,3r-'{Y3,3(f^)-Y3,_3(f^)}, (14) 

where the signs in front of three coefficients Ai^, A^^ 
and A3 3 depend on the sites: the plus signs in front of 
Ai and A3 are taken for A, B, E and F sites and the 
minus signs for the other sites and -1-^.3^3 for B, C, E and 
H sites and —^3^3 for the others. 

The relative values of the these three parameters were 
estimated using a VOg cluster model and they are A^^s/A^fi 
0.79 and {A3fl{4p\r^\3d))/{Ai^o{4p\r\3d)) ~ -0.22. If we ' 
take into account this hybridization as a second order 
perturbation with respect to V'' , we find for the interfer- 
ence term /dq(cli) as 



/DQ(t^) 



q'yRi 



{g\DW)i- 



En 



{?,\Q{e'.q') 



Lo - i7,„ + i^ 



Q{e,q)\g) 



En 



'A 



4p-3d 



)^(^)|g> 



(15) 



where |g) and Eg are the ground state wavefunction 
and its energy, D{e) and Q{e,q) are the operators for 
the electric dipole and quadrupole transitions defined by 
D{£) = e{e-r) and Q{£,q) e{r-£){r-q), A^p-sd is the 
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energy difference between the 4p and 3d levels and iJm 
and Tin are the Hamiltonian and the life-time term for 
the intermediate state reached by the Is — > Sd excitation 
from the initial state. 



The term l/{Eg-^ 
mately replaced by 



Lu — + irni/2) would be approxi- 

'™Y.^i^i^g + - ^^^)\^^){^^\^ where 
^iiiIm) — s-iid with this approximation, the RXS 

amplitude integrated by lo over Is 3d excitation en- 
ergy region can be written as 



Tr{eq) 



2as|r|4p)(ls|r2|3d) 



l.rn,'m' cr 

X [D*^{e')Vi^^„,Q,r.'{e, q) - Q:^{e' , q')Kn',™An'(£) 

(16) 

where V^^^, = J Rip{r)Y*^{n)V\r)R:i^{r)Y,^, (n) dr 
and a|j^^ denotes creation operator of an electron on the 
3c? orbital with magnetic m and spin a quantum numbers 
on site I. Dm and Qm are coefficients of the spherical 
expansions of the transition operators 

1 

m— — 1 
2 



and they are given by 



D±i = 

Q±2 = 

Q±i = 
Qo = 



/ 47r + iSy 




^/2 ' 

— {T{qzex + qxSz) + Kqy^z + q^^v)}: 

lo 



The procedure of the calculations is following. Eight 
V ion sites (A to H in Fig. 11) in the monoclinic unit cell 
was considered to be four independent V ion pairs. For 
each of them, the ground state wavefunction was calcu- 
lated using the V ion pair model described in § ^ In 
the calculations, molecular field Hm — 0.01 eV along the 
magnetic easy axis was applied so to reproduce experi- 
mentally observed antiferromagnetic order. Then, from 
the ground state wavefunctions, the quantity {g\ aimaO- 
were calculated, and using eq. (16), the scattering am- 
plitude was obtained. 

5.3 Analysis of the (111) Bragg reflection 

Figure 12 shows the theoretical azimuthal angle $ de- 
pendence of the integrated intensities for the (111) Bragg 
reflection in the AFI phase for various values of -Dtrg- 
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Fig. 12. Azimuthal angle dependence of the integrated inten- 
sities for the (111) Bragg reflection in the AFI phase calcu- 
lated with various values of the trigonal crystal field -Dtrg- In- 
tensities for the two porlarization conditions a-a and ct-tt are 
shown. The adopted parameter values are -Dtrg = 0.15 eV (a), 
: 0.28 eV (b), Dtrg = 0.335 eV (c) and Dtrg = 0.5 eV (d). 



D 



trg 



The intensities with the two polarization conditions a-a 
and (T-TT calculated with the parameters I?trg = 0.15 eV 
(a), Arg = 0.28 eV (b), Arg = 0.335 eV (c) and Arg = 
0.5 eV (d) are shown. The scattering plane is parallel to 
the z-axis at $ = 0. The tilting angle of the magnetic 
moment in the ground state with -Dtrg = 0.15 eV, 
0.28 eV, 0.335 eV and 0.5 eV are 6»m = 16°, 51°, 71° and 
84°, respectively. 

It is seen in the figure that the azimuthal angle <i> de- 
pendence of the scattering intensities are strongly influ- 
ence by the direction of the orbital magnetic moment. 
The intensity curve in Fig. 12(c) with 9m = 71°, which is 
the same angle observed in the neutron diffraction expcr- 
imentEP and Fig. 12(b) with 6'm = 51° are most resemble 



to the experimental intensity curve in ref. 20, On the 
other hand, the curve in Fig. 12(b), where 9m = 16° and 
the ground state is almost pure e^ai;e'^e'^ {e'^ e'^ ; e'^ ai) 
state, is very different from experimental one. This in- 
dicates that the experimental results are well explained 
with the e'^ e'^ ; e'" e'" and e'^ai;e^e'^ mixed ground state 
and consistent with the present theory. It is also im- 
portant to mention that if the spin-orbit interaction is 
switched off, the scattering intensity is zero for this re- 
flection, although there exists the molecular ferro-orbital 
order in the ground state. The reflection, therefore, orig- 
inates from not the molecular ferro-orbital order but the 

The scattering intensity is also sensitive to a small x- 
component of the orbital magnetic moment M^^- Possi- 
bility of the presence of such small x-component of the 
magnetic moment was pointed out by Moon.LP M^x can 
be induced by the low symmetry crystal field CA(eJ^, e^). 
In Fig. 13(a), $ dependence of the integrated intensi- 
ties for the (111) refiection obtained with CA(eJ^,eJ^) = 
-0.02 eV (sofid line), eV (dashed fine) and 0.02 eV 



Title of the Article 



15 



(a) 



6 

4 

2 


1.5 

1 

0.5 




; A 


(5-n 


















■J \i/^'k^/ ^ 





-160-120-80 -40 40 80 120 160 
Azimuthal Angle (deg.) 

(C) 




Fig. 13. In (a), azimuthal angle dependence of the integrated in- 
tensities for the (111) reflection in the AFI phase calculated with 
CA(eJI,eJ) = -0.02 eV (solid line), 0.0 eV(dashed line) and 
0.02 eV (dotted-dashed line) are shown. In (b) and (c), the 
same as Fig. 11 (c) and (d), the sign (±) of M^^ on V ions are 
shown. 



(dotted-dashed line) are drawn. Z^trg = 0.28 eV is adopted. 
In Figs. 13(b) and (c), the sign of Ml^ on each V ion 
for CA(eJ^,e^) < is shown for the two magnetic lay- 
ers in Figs. 11(b) and (c). For CA(eJ^,e^) > 0, aU the 
signs of Mj^x are altered. As seen in the Fig. 13(a), 
the a-a scattering intensity is strongly influenced by the 
value of CA(eJJ, e^), although the magnitude of Mlx in- 
duced by CA(e^, ej) = ±0.02 eV is an order of- O.OlAiB- 
With CA(e^, ej) — 0.02 eV the calculated intensity curve 
is not resemble to the experimental one whereas with 
CA(eJ^,e^) = —0.02 eV, good agreement with experi- 
mental intensity curve is achieved. 

Lovesey and Knight argued that the (111) reflection 
peak-can be explain by pure Is — s- Sd quadrupole transi- 
tionrS' Contrary to this, in the present theory, the peak 
is attributed to the interference process of the Is Ap 
dipole and Is — + 3(i quadrupole transitions. For compar- 
ison, in Fig. 14, $ dependence of the scattering inten- 
sity of the (111) reflection peak obtained with the pure 
Is — > 3d quadrupole transition instead of the interfer- 
ence process is shown. The intensity curve calculated 
with £)trg = 0.28 eV is very different from experimental 
curve. Results obtained with other parameter values of 
i^trg and CA(e^, e^) also does not agree with experimen- 
tal intensity curve. 
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Fig. 14. Azimuthal angle dependence of the integrated intensity 
for the (111) reflection in the AFI phase calculated with Dtrg = 
0.28 eV. The Is 3d quadrupole transition instead of the Is — > 
3d and Is — > 4p interference process is assumed as for the RXS 
process. 



§6. Linear Dichroism in V 2p XAS 

In this section, the dependence of the V 2p XAS spec- 
tra on the polarization vector e of the light will be dis- 
cussed. As was mentioned in the previous sections, such 
experiment was already done by Park et alr3 with the 
polarization vectors e being parallel and perpendicular to 
the c/i-axis. Because of the selection rule of the 2p 3d 
dipole transition, the transition probability to the oi 
orbital is larger than that to the orbital for ef/ch, 
whereas that to the e'^ orbital is larger for the e _L c/j. 
Thus the difference spectrum is sensitive to the e^ai to 
e'^e'" occupation ratio in the initial state. The ratio was 
deduced from the experiment using a VOg cluster model 
and was found to be e'^ai: e^e^ ~ 1 : 1 for the PM phase 
and e"ai : e"e^ ~ 1 : 2 for the AFI phase. 

To verify the present theory, the spectra was calculated 
using the V ion pair model described in ! 



2.1 



For the ini- 
tial state, the ground state of the V ion pair model is as- 
sumed. For the final state, the configurations 2p'id^\ id"^ , 
2p?>d^;'id^ and 2p3(i^;3(i^ are considered, where 2p de- 
notes a hole on the V 2p orbital in one of the V ion site. 
In addition to the interactions considered in the Hamil- 
tonian described in § 2.1, the 3d-2p multipole interaction 
and the 2p spin-orbit interactioo. are also included in the 
calculation of the final statesO The 2p 3d dipole 
transition is considered as the photoprocess and the in- 
tensities of the spectra (w) as a function of the incident 
photon energy lo is written as 

Uu;)^T\{f\T,\g)\'5{Eg+cj-Ef) 



f 

lim ■ 

(5^0 



Eg + uj - Hf + iS 



where |/), Ef and Hf denote the wavefunction, eigen 
energy and Hamiltonian of the final state, respectively, 
represents the 2p — > 3d dipole transition operator 
and K specifies the polarization direction of the linearly 
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polarized light. Xhe spectra was calculated using the 
recursion methodcj) and a Lorentzian broadening was 
further applied the spectra. 

The parameters used in the calculations are the same 
to those listed in Table II and Dtrg = 0.28 eV is adopted. 
In the final state, there are additional parameters re- 
lated to the 2p core-hole: the 3d-2p Coulomb repul- 
sion energy Udp = 3.07 eV, the 2p spin-orbit interaction 
Cp = 4.649 eV, the 3d-2p Slater integrals F^^ = 4.361 eV, 



= 3.162 eV and G% 



1.797 eV. 
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Fig. 15. Theoretical polarization-dependent V 2p XAS spectra 
for V2O3 for the PM (a) and AFI (b) phases. The spectra are 
depicted for two different directions of the polarization vectors 
e of the light: parallel and perpendicular to the c^-axis. The 
difference spectra are also shown below in each panel. 



Figure 15 shows the calculated spectra obtained using 
the V ion pair model with the polarization vectors e par- 
allel and perpendicular to the c/j-axis. In Fig. 15(a) the 
spectra for the PM phase, and in Fig. 15(b) those for 
the AFI phase are shown and in each panel, the differ- 
ence spectrum is also depicted. The theoretical spectra 
well reproduce experimental ones especially those for the 
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Fig. 16. Theoretical polarization-dependent V 2p XAS spectra 
for the AFI phase calculated without (a) and with (b) the 3d 
spin-orbit interaction. -Dtrg = 0.28 eV is assumed. The spectra 
are depicted for two different directions of the polarization vec- 
tors e of the light: parallel and perpendicular to the monoclinic 
bm, axis. The difference spectra are also shown below in each 
panel. 



AFI phase in ref. |TJ. Although the difference spectrum 
for the PM phase is less similar to the experimental one, 
its tendency to have negative peaks at -7 eV and eV is 
consistent with the experiment. The polarization of the 
spectra defined by P = (/|| -/_l)/(/|| +/_l) are P = 0.027 
for the AFI phase and P = 0.015 for the PM phase, 
where /|| and I± denote integrated intensity of the spec- 
tra with the light polarization parallel and perpendicular 
to the c/i-axis, respectively. Such large change in the po- 
larization between the two phases is also consistent with 
the experiment. Note that although the electron occupa- 
tion of the e"' orbitals in the ground state is small ~ 5%, 
the inclusion of the e" orbital is important to reproduce 
the spectra, especially its polarization. 

The sensitivity of the V 2p XAS spectra with polar- 
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ized light to the orbital occupation in the initial state 
is useful to ascertain whether the orbital order exists or 
not in the AFI phase. To demonstrate this, in Fig. 16 
the polarization-dependent V 2p XAS spectra are shown 
for two different direction of e in the plane which is per- 
pendicular to the c^-axis: e _L b,„ and ej jh^. The 
spectra without (a) and with (b) the Sd spin-orbit in- 
teraction is presented; other parameters are the same to 
those adopted for the AFI phase. The difference spectra 
are also shown below in each panel. 
As was discussed in 



e'^fli; e'^e'^ and 



gTTgTT. gTTgTT jg lai'gc aud thus no electron- 



2.3, in the absence of the 3d 



spin-orbit interaction, the ground state has an orbital 
polarization corresponding to the molecular ferro-orbital 
ordering state of Milla's model (Fig. 16(a)). On the 
other hand, in the presence of the spin-orbit interaction, 
the ground state has no orbital ordering and instead of 
this a large orbital magnetic moment is induced in the 
ground state (Fig. 16(b)). The magnitude of the differ- 
ence spectrum in Fig. 16(a) is much larger than that in 
Fig. 16(b) and the shape of the spectrum in Fig. 16(a) 
largely vary from that in Fig. 16(b). Thus the measure- 
ment of this linear dichroic V 2p spectra, will be a crucial 
check whether the large orbital moment is present or not. 

§7. Conclusions 

The 3d electronic structure and phase transition in 
pure and Cr doped V2O3 were theoretically investigated 
in relation to the id spin-orbit interaction and lattice 
distortion. The 3d electronic state of V2O3 have been 
discussed, using a V ion pair model on the basis of the 
configuration interaction approach. In the model, the 3d- 
3d multipole interaction, a trigonal crystal field and the 
3d spin-orbit interaction were considered and electron 
hopping between nearest neighbor V ions was also taken 
into account. In the presence of the 3d spin-orbit inter- 
action, spin and orbital degrees of freedom are strongly 
coupled and this removes the orbital degeneracy of the 
ground state, resulting in no orbital ordering. Instead of 
this, each V ion with S" = 1 spin state has a large orbital 
magnetic moment ~ 0.7/iB hi the AFI phase. These re- 
sults are in contrast to those obtained from the similar V 
ion pair model proposed by Mila et ai, in which aJeKco 
molecular orbital ordering occurs in the AFI phase Jl3't2P 
The large orbital moment is in accordance with the ratio 
Ml/A/s —0.3 deduced from the nonresonant magnetic 
x-ray scattering experiment)23) and also explains small to- 
tal magnetic moment 1.2/ib observed in the AFI phase. eP 
Because of the presence of the large orbital magnetic mo- 
ment, the monoclinic lattice distortion and the direction 
of the magnetic easy axis is strongly connected in the AFI 
phase. The direction of magnetic moment was found to 
be canted from the corundum c^-axis in the presence of 
the monoclinic distortion and this agrees with the neu- 
tron diffraction experiments.^ 

The local lattice distortion energy due to tilting of the 
V ion pair from the corundum c-axis was further consid- 
ered. From the present model, it is expected that there 
is at least two local minima for the free energy in the 
real system as a function of the V-V distance of the pair. 
One is positioned at the V-V distance in the PM phase, 
where the energy difference of the two kinds of the state 



lattice coupling takes place. The electronic state of the 
V ion pair in this phase is described as a superposition 
of two configurations e'^ai]e'^e^ and e'^e'^;e'^ai in equal 
weight, which is similar to Mila's model. The other one 
is located at the V-V distance in the PI and AFI phases, 
where energies of the two kinds of the states are degen- 
erate and Jahn- Teller like lattice distortion is the cause 
of the energy lowering at this V-V distance. Because of 
the tilting of the V ion pair caused by this lattice in- 
stability, the e^e^;e^e^ state is further hybridized with 
above e^ai; e^e^ (e^e^; e^ai) state in these phases. This 
large difference in the electron occupation between the 
metal and insulating phases agrees with the recent lin- 
ear dichroic V 2p x-ra^^ absorption spectroscopy (XAS) 
experimental results.ES 

Since both in the AFI and PI phases, the electron- 
lattice coupling strongly influences the orbital occupa- 
tion and charge fluctuation, the PM AFI and PM 
PI transitions can not be regarded as usual Mott transi- 
tion. The elongation of the V-V distance and resultant 
reduction in the hybridization strength between the 3d 
orbitals caused by this electron-lattice coupling is prob- 
ably responsible for these metal-insulator transitions. 

To investigate interplay between magnetic ordering and 
the lattice distortion, an effective spin Hamiltonian in- 
cluding the effects of tilting of the V ion pairs was intro- 
duced. The antiferromagnetic order with the monoclinic 
lattice distortion observed in the AFI phase was repro- 
duced in this model, where simultaneous ordering of the 
tilting directions and the magnetic moments of the V 
ion pairs take place. With increasing temperature, it 
exhibits antiferromagnetic to paramagnetic phase tran- 
sition corresponding to the AFI PI transition in the 
real system. While all V ion pairs are tilted to the same 
direction and this gives rise the monoclinic lattice dis- 
tortion in the AFI phase, their tilting directions are dis- 
ordered and fluctuate among the three stable directions 
in the PI phase. As a result, there is no monoclinic dis- 
tortion in the PI phase. It was also found that the C44 
elastic constant softens and short-range spin correlation 
functions abruptly changes at this second order transi- 
tion in this model. These resuLts, are consistent with 
the ultrasonic-wave meiasurementtlP and recent neutron 
scattering experimentsBll^l* 

To confirm the validity of the present model, recent 
experiments OH_the resonant X-ray scattering (RXS) at 
the V X-ed^BEJ and the hnear dichroic V 2p XAS 
experimentsEJ' were analyzed. The (111) Bragg reflec- 
tion observed in the RXS experiments, particularly in 
its azimuthal angle and polarization dependence of the 
Is — > 3d peak, is well explained within the present model, 
where no orbital ordering is present and the antiferro- 
magnetic order is assumed. The reflection is pure mag- 
netic and the scattering amplitude of this reflection mainly 
arises from the interference process of the Is ^ 4p dipole 
and Is —> 3d quadrupole transitions. The shape and 
magnitude of the linear dichroic V 2p XAS spectra in 
both the AFI and PM phases are well reproduced within 
the present model. Another linear dichroic V 2p XAS 
experiment has been proposed, which is sensitive to the 
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presence of the orbital order in the initial state and will 
be a crucial check to the present model. 
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